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| Abstract 

| Known perturbative expressions for the decay rates of 115 GeV Higgs 

D . particle into either two gluons or a bb pair are shown to exhibit renormalization- 

scale- (^)-dependence that is largely removed via renormalization-group/Pade- 
approximant estimates of these rates' next order contributions. The ex- 
trema in fi characterizing both rates, as determined from fully-known 
orders of perturbation theory, are very nearly equal to corresponding /i- 
insensitive rates obtained via estimation of their next order contributions, 
consistent with "minimal-sensitivity" expectations. 

Physical hadronic processes, such as measurable decay rates, are necessarily 
independent of the renormalization-scale parameter [i that characterises order- 
£\j ' by-order calculations within perturbative QCD. Significant renormalization- 

scale dependence, however, does appear within QCD calculations that are trun- 
CD ' cated after only two or three orders of perturbation theory. Since the Feynman- 

diagrammatic calculation of subsequent orders of QCD is very often prohibitively 
' difficult, a common practise for decay rates is to identify /i with the mass of 

Oh, the decaying particle. Alternatively, one might either choose a value of /i for 

O 1 which the /z-dependence of the known terms of the rate is locally flat, reflecting 

D , (hopefully) the renormalization-scale invariance of the all-orders calculation, or 

one might choose a value of /i for which the highest-order known term of the 
perturbative series is minimized, thereby forcing an (apparent) convergence of 
the known terms of the perturbative series. 

All such procedures for identifying [i, however, are necessarily prescriptions 
for obtaining a prediction, and any discrepancies between predictions from dif- 
fering procedures necessarily contributes to the theoretical uncertainty of such 
predictions. One purpose of the present work is to demonstrate how this (of- 
ten overlooked) component of the theoretical uncertainty can be reduced via 
estimation of the next-order of the perturbative series from which theoretical 
predictions are obtained. 

In the present work, we consider the two gluon (gg) and bb decay modes 
of a Standard Model Higgs particle (H), as obtained from spontaneous elec- 
troweak symmetry breaking via a single complex scalar- field doublet. We as- 
sume that Mh = 115 GeV, as suggested by recent exciting results emerging from 
LEP Q. The non-mass-suppressed contributions to the H — > gg and H — ► bb 
rates have been respectively calculated to two and three nonleading orders of 
perturbative QCD. We show here that these calculated rates indeed exhibit 
renormalization-scale dependence, although this dependence is more pronounced 
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for H — ► gg than for H — ► £>&. We discuss how the next term in the pertur- 
bative series for both processes can be estimated via asymptotic Pade approx- 
imant methods developed by Samuel and collaborators |^|, ||. We demonstrate 
the success of such methods in predicting those next-order terms accessible via 
renormalization-group methods, and show how the incorporation of the esti- 
mated next-order term dramatically reduces the theoretical uncertainty associ- 
ated with renormalization-scale dependence for both processes. We then obtain 
predictions for H — > gg and H — ► bb decay rates for a 115 GeV Higgs particle 
for which those theoretical uncertainties associated with renormalization-scale 
dependence are virtually eliminated. 

The Higgs — > gg decay rate, as extracted from correlation functions in the 
large (2M t /M H ) 2 limit (§, may be expressed in the following form Q 



H^gg 



V2G F Mfj 



72tt 



S(jj) 



(1) 



where the renormalization-scale-sensitive factor in the decay rate (1) is given by 



S(n) = S[x(n),L(n),T\ 
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with 



xOji) = aPM/w, (3) 

L(jm) ee \og(v 2 /m 2 t (ri), (4) 

T = \og(M 2 H /M 2 ) pole . (5) 

Equation (2) reflects the fact that the decay rate is known to three-loop order, 
and includes an appropriate parameterization of the subsequent unknown four- 
loop-order contribution. 

For the decay rate of a 115 GeV Higgs particle into two gluons, this scale- 
sensitive factor is completely determined at /i = Mh by the known central 
values of ai 6 \M z ) = 0.1185 GeV (from which coupling-constant evolution can 
proceed to arbitarary fi) and M t — 174.3 GeV Q. Expressed in terms of its 
three known leading orders, this scale-sensitive factor (2) is given by 



S(M H )3-ioop = 0.00248 = (0.1148/tt) 2 [1 + 0.657 + 0.201] 



(6) 



Because of its relatively large final term, the truncated perturbative series within 
(6) appears likely to underestimate the true value of the series sum. 



2 



The factor S(/j.) satisfies the renormalization group (RG) equation 



dS 



dp 2 



= p 2 —[ x ( fi ),L(p),T} = l-2^(x) — +^(x) — , (7) 



dS ,dS 



dL dx 



a reflection of the renormalization-scale invariance of the physical decay rate. 
Perturbative expansions of the right-hand side of (7) can be utilized to obtain 
order- by-order cancellation of scale-sensitivity within S(/i). Thus the value of y, 
chosen in (6) should, in principle, not affect the rate obtained via (1). However, 
truncation of the perturbation series within (2) and (6) necessarily introduces 
a residual scale-dependent within the rate. In Figure 1, we have plotted the 
three-loop expression for S(p) for /i over the entire m&(m&) < fi < M t range 
of physical interest. The three loop rate is manifestly renormalization-scale 
sensitive at /i = 115 GeV, exhibiting an increase with decreasing \i until local 
scale-insensitivity is obtained near /i = 49.5 GeV: 

S(49.5 GeV) 3 - Loop = 0.00254 = (0.1293/tt) 2 [1 + 0.475 + 0.026]. (8) 

Note that the convergence of the series within (8) is much more evident than 
within (6), reinforcing the suggestion that (6) underestimates the full series sum 
at \x = Mjj- 

In previous work [^J, we have demonstrated how asymptotic Pade approx- 
imant methods can predict the four- loop coefficients Cj(T) within (2), where 
T is explicitly given by the renormalization-scale invariant logarithm (5): T = 
21og(115.0/174.3) = -0.8317. Thus, the perturbative series within (2) is of the 
form 

S(p)/x 2 (p) = 1 + Ri(L)x + R 2 (L)x 2 + R 3 (L)x 3 + (9) 

where 

Rx{L) =21.10+ -L (10) 
147 

R 2 (L) = 239.5 + 107.LL+ L (11) 

16 

Ra(L) = c + Cl L + c 2 L 2 + c 3 L 3 (12) 
The asymptotic Pade-approximant prediction , 

Rf + R1R2 



^3 * ,,3 , I n ■ (13) 



as derived M from the anticipated error g of [N\l] approximants in predicting 
the term Rn + 2X N+2 within the series (9), can be matched to the form of (12) 
by equating the four lowest moments of both versions of R3 : 

w = m 2 t (fx)/fx 2 , (14) 



(k + 2) / dw w k+1 [c — ci log w + c 2 log 2 w — C3 log 3 w] 
Jo 



(k + 2) / dww* 



k+l 



(i 



2R*(L) 



R{(L) + R 1 (L)R 2 (L) 



, (L = -log(w)).(W) 
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The final line of (15) is evaluated numerically via explicit substitution of (10) 
and (11) into the integrand of (15), leading to the following four equations: 



with solutions 



N-i = c + ci + 2c 2 + 6c 3 = 4267 

N = c + a/2 + c 2 /2 + 3c 3 /4 = 2858 
Ni = co + ci/3 + 2c 2 /9 + 2c 3 /9 = 2493 
N 2 = co + d/i + c 2 /8 + 3c 3 /32 = 2329 

c = 1900, a = 1527, c 2 = 359.3, c 3 = 20.20 



(16) 

(17) 
(18) 
(19) 

(20) 



The accuracy of these solutions is, in part, testable by comparison of the 
predicted values for {ci,c 2 ,c 3 } with the values one can extract from the RG 
equation (7) . Given the perturbative expansions of nf =6 RG functions 



/3 (6) (x) = -x 2 
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x 2 + 0{x 3 ) 



(21) 
(22) 



the coefficients c±, c 2 , and c 3 are easily obtained by the requirement that the 
aggregate coefficients of a; 5 , x 5 L, and x 5 L 2 on the right-hand side of the RG 
equation (7) all vanish: 



„RG 



= 1540.4, cf G = 364.83 



cf G = 21.4375 



(23) 



The striking agreement between the true values (23) and the estimated values 
(20) suggests that the estimate in (20) of the RG-inaccessible coefficient c is 
comparably accurate. One test of the internal consistency and stability of the 
estimation procedure is to substitute the known values (23) for {ci, c 2 , c 3 } into 
eqs. (16-19), and then to solve each of these equations for q> 

iV_i : co = 1868 (24) 

(25) 
(26) 
(27) 

All four of these values are in good agreement with the estimate Co = 1900 
[eq. (20)] obtained without the use of RG values for {ci, c 2 , c 3 }. Moreover, the 
estimate (24) also follows from optimization dy? /dco = of the least-squares 
fit of (12) to (13) over the UV domain of /j 2 , as discussed in S: 
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* G log(w) 
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Thus, to estimate the four loop contribution to we choose to incor- 

porate the estimate (24) for cq in conjunction with the known values (23) for 
{ci,C2,cz}. We also incorporate leading mass corrections to the rate, which 
can be expressed as a further contribution 5 m cc 2 (/i) to S([i), where || 

S m = ^ [(A t + ReA b ) 2 + (ImA b ) 2 ] - 1, (29) 

A t =2[r t + {n - l){sin-\^)) 2 } /r 2 , T t = M 2 H /AM 2 (30) 



A b = 2[T b + (T b -l)f(T b )]/T 2 , 

U = M 2 /4m 2 , 



f(r) = -\ 



lOff L ^=^= ITT 

8 (1-0^177) 



(31) 



This correction to S(fi) is quite small compared to the contributions of the 
leading three-orders of perturbation theory, but is comparable in size to the 
four loop term. For example, using the "physical" mass values M t — 174.3 
GeV, M H = 115 GcV, and m b = 4.20 GeV we find that S m = -0.0597, and 
that the mass-corrected four-loop version of (6) is given by 

S(M H ) 4 - loop+s = 0.00245 = (0.1148/7r) 2 [l + 0.657 + 0.201 + 0.0373 (-0.0597)]. 

(32) 

In (32), the final number in parantheses is the mass-correction factor 5 m , pre- 
ceded by the one-loop (unity), two- loop, three-loop and (estimated) four- loop 
contributions to the rate in the AM 2 /M^j >> 1 limit, as obtained via (2), 
(23) and (24). These series terms exhibit on accelerated convergence over the 
four-orders of the perturbation series within (32); the ratio of successive terms 
(respectively 66%, 31% and 18.5%) is itself seen to decrease as the order in- 
creases. More importantly, the mass-corrected four-loop rate 5*(/i) is virtually 
flat over a very broad range of /i inclusive of fi = Mh, as is evident from the 
"4L" curve displayed in Figure 2. This curve is seen to exhibit a very weak 
local maximum at [i = 110.5 GeV which differs infinitesimally from the S(Mh) 
rate of (32); i.e. [5(110.5 GeV)/S(ll5 GeV)] A _ loop+Sm = 1.00002. Had we 
not included the mass correction factor 5 m , this point of minimal sensitivity to 
renormalization-scale would have occurred at 74.5 GeV. 

In Figure 2, we have also displayed the mass-corrected version of the (fully- 
known) three-loop (3L) expression for S(fi); for both the 3L and 4L curves in the 
Figure, the leading-order mass correction S m x 2 (^) is included in S(/i). The 3L 
curve clearly exhibits more renormalization-scale dependence than the 4L curve. 
However, the magnitude of the local maximum of the 3L curve (occurring at 
[i = 57.5 GeV) is in very good agreement with the plateau portion of the 4L 
curve, i.e. the estimated four- loop rate over the broad region of insensitivity to 
the renormalization-scale ijl. Such agreement of the 3L curve's extremum with 
the scale-insensitive next-order result is anticipated by arguments presented by 
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Stevenson two decades ago Q, in which the physical result to a given order 
of perturbation theory is identified with the point of local insensitivity to an 
(ultimately) unphysical parameter. The near equivalence of the 3L curve's ex- 
tremum with the 4L curve's plateau supports the prescription of choosing a 
renomalization scale for which a perturbatively-obtained expression is "mini- 
mally sensitive" p0| , i.e. locally flat. Thus, the factor S(n) within the rate 
(1) is estimated to four orders of perturbation theory by (32), corresponding to 
T H ^ gg = 2.71 • 1CT 4 GeV if M H = 115 GeV. As is evident from Fig. 2, this 
estimate is insensitive to the magnitude of the renormalization- scale parame- 
ter over a very broad domain of /i beginning at ji = 20 GeV and inclusive of 
/i = Mh, in contrast to the known three-loop expression for S(/i). 

Within a Standard Model single-Higgs-doublet context, however, the process 
H — > bb is the dominant hadronic decay mode of a Higgs particle with mass 115 
GeV. This decay rate has been obtained explicitly by Chetyrkin |ll| , and may be 
expressed in terms of a renormalization-scale dependent factor i?(/i) as follows: 

T(H -> bb) = ^f^(^) (33) 

where, for a = ai 5 \fj,)/n, the scale-dependent factor R(fi) is proportional to 
the sum of H[a(/i), /i], the imaginary part of the scalar correlation function, and 
Higgs-mass-suppressed contributions £[a(/i), m;,(/i)]: 

R(fi) = ml(p) PI[o(aO, M] + S[o(m), m b (ji)]] , (34) 



Il[a(fi,fi] = 1 + a(fi) 



"!7 „, 






y + 21og 





29.14671 



a 2 ( M ) 

« 3 (M) 
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di log 4 



£[a(/z), m 6 (/i)] = -g^ [-6 - 40 a( M ) - 87.72 a 2 (,i) + ...]. (36) 

The factors of Mh appearing in (33), (35) and (36) are pole masses that are 
independent of the renormalization scale parameter fi. In (35), we have in- 
cluded explicit parameterization of the C(a 4 ) contribution to II. The coefficients 
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{d\, d,2i 1^4} can be extracted via the RG equation 113] 

M 2 ^2 [ m b(M) n [ a (M)>M]] =0 

n[a(/*),d=0, (37) 



so as to yield the following values ]12[ ]: 

di = 791.5233, d 2 = 1114.695, d 3 = 260.0647, d 4 = 14.75550. (38) 

However, the non-logarithmic coefficient do is inaccessible via RG arguments at 
this order; its value can be obtained only by an explicit Feynman diagrammatic 
calculation which, as of the present time, has not been performed. 

Earlier work ^ |7[ ^) has explored asymptotic Pade-approximant approaches 
toward the estimation of this do term. Given a field theoretical series of the 
form 1 + aRi + a 2 i?2 + cl 3 R 3 + a 4 i?4 + where {Ri, R2, R3} are known and 
i?4 is unknown, the 0(l|iV) anticipated error of [iV|l] Pade-approximants in 
reproducing the series coefficient Rn+2 [|| can be utilized to predict that |Q] 

j^pred. _ RljRy + R2R2R3 - 2R 3 R 3 ] 

i? 2 [2i?2 — R\Ri — ^1^1] 
If, as in the pseudoscalar correlator (35), the coefficient of a 4 is given by 

H 4 = d Q + dl log + d 2 log* [-Q + d 3 log3 ( J^) + d 4 log 4 ( 

. (40) 

the coefficients {do, di, d2, d 3 , d^} can be predicted by an appropriate fitting of 
(39) to the form of (40). The results of such fitting jlj] are predicted values 
[di = 745, d 2 = 1180, d 3 = 253, d 4 = 15.4] within ±6% of the RG-determined 
values in (38), as well as a prediction that do = 64.2. This prediction appears 
to be corroborated by the (rif = 5) estimate do = 62 recently obtained by 
Broadhurst, Kataev, and Maxwell [|l3|, as well as an earlier estimate (do = 67.2) 
obtained via (39) entirely from nonlogarithmic terms [Q. 

In Figure 3 we have plotted R{^), the renormalization-scale sensitive portion 
of the H — > bb decay rate for 

1) the case (4L) in which II, as given by (35), is truncated after its fully-known 

0(a 3 ) contribution, and 

2) the case (5L) in which the 0(a A ) contributions to the correlator (35) are ex- 

plicitly incorporated through use of the RG values (38) for the coefficients 
{di, d2, d3, d 4 } in conjunction with the estimate do = 64 obtained in p2[ 
via asymptotic Pade-approximant methods. 

For both cases, the same mass-suppressed contribution £[a(/i), m,b(fj.)] is 
employed. The final 0{a 2 ) term of this contribution, as given in (36), is generally 
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smaller than the estimated 0(a 4 ) contribution to the pseudoscalar correlator 
(35). In obtaining both Figure 3 curves, the evolution of a(fi) from the central 
empirical value || a s {Mz) — 0.1185 is obtained via the known terms of the 
n/ = 5 /3-function Q: 

2 da 23 2 29 3 9769 4 1Q oroi „ 5 
"V = ^ - 12° ~ 3456° ~ 18 - 85217a " (41) 

Similarly, the evolution of rrib(fi) from the assumed five-flavour threshold at 
mf,(4.2 GeV) — 4.2 GeV is governed for both curves by the known terms of the 
rif = 5 7 m -function Jl5| : 

/ \ .o^t/T a ^ 1 12/23 C K^)] 
mb(/i) = 4.2 GeV 



a(4.2 GeV) 



c[a(4.2 GeV)} ' 



c[a] = 1 + 1.17549a + 1.50071a 2 - 0.172486a 3 . (42) 

A comparison of the two curves shows that inclusion of the estimated C(a 4 ) 
contribution to II leads to a significant flattening in the renormalization-scale 
dependence of the calculated H — > bb rate. Once again, the extremum of i?4i(/x), 
the curve for which estimated 0(a 4 ) contributions to II are excluded, corresponds 
to the broad 9.85 GeV 2 plateau in i?5^(/i), the curve which does incorporate 
the estimated C(a 4 ) contributions to II, consistent with the expectation that 
the point of minimal sensitivity to /i should correspond most closely to the true 
rate jlO| . The extent of this plateau in R^l (/i) can be ascertained by noting 
that R$l{h) is between 9.840 and 9.860 GeV 2 over a very broad domain of n 
(22.5 GeV < /z < 313 GeV). By contrast, the values of i?4i(/x) are within the 
same range only for a substantially smaller domain of /i (69.5 GeV < fi < 147 
GeV). Thus, it is evident that the inclusion of the asymptotic-Pade-approximant 
estimate of the next order contribution to the scalar correlator significantly 
reduces the renormalization-scale dependence of the predicted H — > bb decay 
rate. If we identify the true rate with the explicit plateau-region value i?5i(/x = 
115 GeV) = 9.851 GeV 2 , the corresponding (M H = 115 GeV) value for the 
decay rate is T(H bb) = 2.231 ■ 10~ 3 GeV. 

We conclude by noting for both H — > gg and H — > bb decay processes 
that if one considers only the fully-known terms of the appropriate perturbative 
series, our findings support the "minimal sensitivity" prescription [1^| in which 
/i is assigned that value for which the rate is locally flat. The decay rates 
obtained from fully-known terms by this procedure are found to be very nearly 
the same as those scale-insensitive rates we obtain via estimation of the next- 
order contribution to the perturbative series. 
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Figure 1: The scale-sensitive portion [5(/x)] of the three- loop order Higgs — > two 
gluon decay rate, as given by (1) and (2), is plotted for M t = 174.3 GeV and 
Mh = 115 GeV, with the scale-dependent coupling evolving from a s (Mz) = 
0.1185. The renormalization-scale parameter \i is in GeV; S(fi) is dimensionless. 
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Figure 2: The scale-sensitive portion S(fi) of the Higgs — > two gluon decay rate, 
as obtained from fully- known contributions up to three- loop-order (3L), is com- 
pared to S(fi) as obtained from inclusion of its estimated next-order (4L) con- 
tributions. To isolate the effect of incorporating the next-order estimate, both 
expressions for S(/j) are inclusive of the mass correction <5 m x 2 (/z), as obtained 
from (29-31) in the text, with m b = 4.2 GeV, M t = 174.3 GeV, Mjj = 115 GeV, 
and the scale-dependent coupling evolving from a s (Mz) = 0.1185. The 
renormalization-scale parameter \x is in GeV; S(fi) is dimensionless. 
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Figure 3: The scale-sensitive portion R(p) of the Higgs — > bb decay rate, as ob- 
tained from its fully-known four-loop-order (4L) contributions to the correlation 
function (35), is compared to as obtained from inclusion of the scalar cor- 
relator's estimated next-order (5L) contributions. As before, M t = 174.3 GeV, 
Mh = 115 GeV. The scale-dependent coupling and b-quark mass are evolved 
from a s (Mz) — 0.1185, rrib(mb) = 4.2 GeV. The renormalization-scale param- 
eter [i is in GeV, and the vertical axis is in GeV 2 units. 
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